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We derive the quantum Boltzmann equation for the two-dimensional electron gas in a magnetic 
field such that the filling factor v ^ 1. This equation describes all of the effects of the external 
fields on the impurity collision integral including Shubnikov-de Ifaas oscillations, smooth part of the 
magnetoresistance, and non-linear transport. Furthermore, we obtain quantitative results for the 
effect of the external microwave radiation on the linear and non-linear dc transport in the system. 
Our findings are relevant for the description of the oscillating resistivity discovered by Zudov et ai, 
zero-resistance state discovered by Mani et al. and Zudov et at, and for the microscopic justification 
of the model of Andreev et al.. We also present semiclassical picture for the qualitative consideration 
of the effects of the applied field on the collision integral. 
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I. INTRODUCTION 

The purpose of this paper is to construct the theory 
describing linear magnetotransport, non-linear effect of 
dc electric field, and effect of microwave on both linear 
and nonlinear dc magnetotransport from a unified point 
of view. Despite a long history of the systematic exper- 
imental and theoretical study of the properties of two- 
dimensional electron and hole systems in semiclassically 
strong-'^ and quantizing magnetic field^, the system still 
brings us surprises. 

Recent experiments'^"^ revealed the new class of phe- 
nomena. (In fact, such effects were first considered the- 
oretically by Ryzhii^'^° three decades ago but were not 
fully appreciated.) Exposing the two-dimensional elec- 
tron system to a microwave radiation, Zudov et al.^ dis- 
covered the drastic oscillations of the longitudinal resis- 
tivity as a function of the magnetic field. The period 
of these oscillations was controlled only by the ratio of 
the microwave frequency to to the cyclotron frequency 
bJc. Moreover, the oscillations were observed at relatively 
high temperature, T, such that usual Shubnikov-de Haas 
oscillations in the absence of the microwave irradiation 
were not seen, T ^ fiojc. Working with cleaner sam- 
ples, almost simultaneously, two experimental groups*'^ 
reported observations of a novel zero-resistance state in 
two-dimensional electron systems, appearing when the 
oscillations of the resistivity hit zero. It is worth em- 
phasizing that the zero-resistance state was not con- 
nected to any significant features in the Hall resistiv- 
ity in contrast to that for the Quantum Hall Effect^. 
Further experimental activity consisted in analysis of 
the low-field part of the oscillations in order to under- 
stand the effect of the spin-orbit interaction^ and obser- 
vation of the zero-conductance state in the Corbino disk 
geometry^. Results^'^ were later confirmed by an inde- 
pendent experiment.® 

Two recent theoretical papers^^'^^ are likely to explain 



the main qualitative features of the data'^"'*. Durst et 
al.^^ presented a physical picture and a calculation of 
the effect of microwave radiation on the impurity scat- 
tering processes of a two dimensional electron gas [see 
also Ref. 10]. In addition to obtaining big oscillations of 
the magnetoresistance with the right period of Ref. 3. the 
crucial result of Refs. 9-11 is the existence of the regimes 
of magnetic field and applied microwave power for which 
the longitudinal linear response conductivity is negative, 

(Jxx < 0. (1.1) 

It was shown by Andreev et al.^^ that Eq. (1.1) by 
itself suffices to explain the zero-dc-resistance state ob- 
served in Refs. 4,5, independent of the details of the mi- 
croscopic mechanism which gives rise to Eq. (1.1). The 
essence of Andreev et al.^^ result is that a negative lin- 
ear response conductance implies that the zero current 
state is intrinsically unstable: the system spontaneously 
develops a non-vanishing local current density, which al- 
most everywhere has a specific magnitude jo determined 
by the condition that the component of electric field par- 
allel to the local current vanishes, see also Sec. VII of 
the present paper. The existence of this instability was 
shown to be the origin of the observed zero resistance 
state. It is worth mentioning that the instability of the 
system^^ with absolute negative conductivity is known 
since the work of Zakharov.^^ The important new feature 
of the instability and the domain structure of Andreev et 
al.^^, that the instability occurs at large Hall angle; as 
the result the domains for the current coincide with the 
domains of the electric field directed perpendicular to the 
current. We also would like to notice the certain simi- 
larity with the model of photoinduced domains proposed 
by D'yakonov^*^ as an explanation of the experiments on 
ruby crystals under the intense laser irradiation 

Subsequent theoretical works outlined ideas^^"^^ of 
Ref. 11,12; postulated^° the plasma drift instability; con- 
sidered "a simple classical model for the negative dc con- 
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ductivity" due to non-parabolicity of the spectmm^"'^ or 
due to the lattice effects on ac-driven 2D electrons. We 
will not discuss those works further in a present paper. 

Unfortunately, the comprehensive quantitative de- 
scription of the data''^"'^ is not possible within Refs. 9- 
11. Moreover, phenomenology of Refs. 12 implies certain 
form of the non-linear dc transport in the presence of 
the microwave radiation which has not been microscop- 
ically justified yet. Our article presents a program for 
such a description. However, we will not take into ac- 
count effects which depend on the distribution function 
and are determined by inelastic processes, see Ref.^^ and 
Section II A. These effects will be considered elsewhere. 

The paper is organized as follows. Qualitative discus- 
sion based on a consideration of semiclassical periodic 
orbits is presented in Sec. II. The quantum Boltzmann 
equation applicable for large filling factors and small an- 
gle scattering on the impurity potential is derived in 
Sec. III. This equation is later used to obtain closed 
analytic formulas for linear-fio transport. Sec. IV, non- 
linear dc transport, Sec. V, and the effect of microwave 
radiation on the dc transport. Sec. VI. Section VII re- 
lates the results to the model of domains of Ref. 12. Our 
findings are summarized in Conclusions. 



II. QUALITATIVE DISCUSSION 

The qualitative discussion of the effect of the mi- 
crowave radiation on the dc-transport was presented in 
Refs. 9-11 in terms of quantum transitions between Lan- 
dau levels. We chose to utilize the fact that only electrons 
with large Landau level indices are important and ex- 
plain the effects in terms of semiclassical periodic motion. 
This explanation becomes especially convenient when the 
Landau levels are significantly broadened, which means 
that the number of the repetitions in the periodic orbit 
is small. [Infinite number of the repetitions of the pe- 
riodic orbit would correspond to the vanishing width of 
the Landau levels.] Moreover, the qualitative picture will 
enable us to separate effects into two groups according to 
their sensitivity to the electron distribution function, and 
understand the status and validity of the approximation 
which will be made in the technical part of the paper. 

To analyze the effect of external fields on the collision 
processes, it is more convenient to switch into the moving 
coordinate frame 

r^r-CW, (2.1) 

in which the external electric field is absent. The position 
of the moving frame <^{t) is found from 

where E{t) is the applied spatially homogeneous electric 
field, me is the electron band mass, ajn = is the 



cyclotron frequency, and e is the antisymmetric tensor: 

^xy — ^yx — 1-1 ^xx — ^yy — 

If there were no disorder potential, the distribution 
function /(e) of the electrons in this moving frame would 
be the Fermi function 

/(^) = M^) = i:p^' (2-3) 

no excitations would appear and therefore no dissipative 
current would be possible. On the classical level, an elec- 
tron experiences the cyclotron motion with the position 
of the guiding center R intact. Collision with impurities 
moving with velocity dtC, cause the drift of the guiding 
center, so that the current density 

arises. Here, A^e is the electron density, and {dR/dt)^^^ 
symbolizes the probabilistic change in the position of the 
guiding center to be discussed below. 

Let us consider the scattering process of the electron off 
one impurity. Because the size of the scattering region 
(correlation length of the potential ^) is much smaller 
than the cyclotron radius, we can still characterize the 
scattering process by the initial direction i, and the scat- 
tering angle 9 as shown on Fig. 1. We consider only small 
angle scattering 

9^1. (2.5) 

In this case, each scattering event causes the shift in the 
position of the guiding center 

AR = -iR^9, (2.6) 

where Rc = vf/oJc is the cyclotron radius and vp is the 
Fermi velocity, see Fig. 1. 

During the collision process impurity moves with the 
velocity —dtC{t)- Because the size of the impurity is 
small, and we assume that |9(^(t)| <C vp, this motion 
can be neglected in the calculation of the scattering am- 
plitude but it has to be taken into account in the conser- 
vation of energy. Indeed, during the scattering event the 
moving impurity transfers the energy Ae = —dtC{t) ■ Ap, 
where change of the electron momentum is given by 

Ap = -sippO (2.7) 

for 9 <ti 1, and pf = meVp is the Fermi momentum. 

Taking this energy change into account we write for 
the displacement of the center of orbit 

= {j d9^ j de{f{e)-f{e + Ae))M). (2.8) 
= lijd99j ds^^ {dtm ■ Ap)m)., 



FIG. 1: Scattering process off a singie impurity. Inset sliows 
tlie semiclassical trajectory in the vicinity of the impurity. 



where the function [9) is proportional to the scatter- 
ing cross-section and is determined by the impurity po- 
tential. We will assume that M.{9) vanishes rapidly at 
9 - h/{pFC) < 1, i-e- Eq. (2.5) holds. In Eq. (2.8), 
{■ ■ ■)l stands for the averaging over the direction of the 
momentum of the electron incoming on the impurity. 

Substituting Eq. (2.7) into Eq. (2.8), one finds 



where 



1 _ Pfvf 

Ttr 2 



VF 



(2.9) 



d99^M{9) 



(2.10) 



is the transport scattering time at zero magnetic field. 
Together with Eq. (2.4) this gives the Drude formula 



for the large Hall angle uJcTtr 3> 1. 



It is not the end of the story though. Considering one 
scattering event as a complete real process, we imply that 
there are no returns of an electron to the same impurity, 
or, to be more precise, the possible returns are not corre- 
lated with original scattering. However, in magnetic field 
an electron moves along a circle of the cyclotron radius 
Rc between shattering processes. This circular motion 
results in correlated returns of the electron to the same 
impurity. [To the best of our knowledge the first dis- 
cussion of the magnetoresistance in terms of returning 
semiclassical orbits was performed in Ref. 25] 

Such returns do not change the structure of Eq. (2.8), 
but they do change the scattering cross-section Ai [9) in 
comparison with its value in zero magnetic field. Indeed, 
one can see from Fig. 2 that several semiclassical paths 
characterized by different number of rotations and dif- 
ferent instances of the impurity scattering contribute to 
the same final state; amplitudes for such processes Af{9) 



FIG. 2: Different amplitude giving coherent contributions into 
the impurity scattering cross-sections. Inset shows the shift 
of the impurity between the scattering events in moving co- 
ordinate frame. 



sum up coherently. 



M{9) cx 



= iAo(0)iV2Re Ar;{9)[Ar{9)r , 

(2.11) 

where index I labels the number of rotation and a labels 
semiclassical paths; Fig. 2 shows the paths for Z = 1,2, 
while Fig. 1 depicts the path for / — 0. Equation (2.8) 
takes into account the contribution from the shortest tra- 
jectory (first term in the second line) of Eq. (2.11) and 
misses the interference contributions. 

To assess the role of the interference contributions, let 
us employ the Born approximation of the impurity scat- 
tering. Then, each semiclassical path may involve only 
one scattering off an impurity, and all the paths are clas- 
sified by (i) the scattering angle 9 and (ii) whether the 
impurity affects the electron in the beginning or in the 
end of the path; we will call corresponding amplitudes 
^(9) and Af[G), see Fig. 2. Factorizing the impurity 
scattering potential into the scattering cross-section at 
zero magnetic field M.q{&), we obtain 

M{9)=Ma{9) \ l + 2ReEyl>(6') {A\ 

(2.12) 

Equation (2.12) neglects the motion of the impurity 
during the whole collision process. Apparently, it is con- 
sistent with the derivation of Eq. (2.8) where the effect 
of the impurity motion on the matrix elements was also 
neglected. Indeed, for the process in Fig. 1, Z = 0, 
the characteristic scattering time can be estimated as 
At — (,/vf- The displacement of the impurity is, then, 
|Ar| = l^tCI • {i/vp) < e, for l^tCI < VF, and could be 
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neglected. However, for the scattering shown in Fig. 2 
the interfering processes are separated in time by the 
interval At = Displacement during such interval 

is A|r| ~ IdtCl ■ (^) = \dtC\ (^) (2^). Because 
Rc^ the displacement may easily become comparable 
with the impurity size ^ even for <C vp- Therefore 
the displacement must be taken into account in the inter- 
ference terms in Eq. (2.12). Due to the impurity motion, 
the scattering off impurity occurs at different points, see 
inset in Fig. 2. In this respect, the scattering off a mov- 
ing impurity is analogues to the interference in a grate 
interferometer with the distance between slits equal to 

Ci{t) = C{t-^)-C (*)• By analogy with the grate in- 
terferometer, we find that Eq. (2.12) has to be modified 

as 



where Ap is given by Eq. (2.7). Finally, the quantum 
mechanical amplitudes Af'^(9) can be decomposed into 
the smooth prefactor a and the rapidly oscillating expo- 
nent involving the reduced action along a semiclassical 
cyclotron trajectory 



af''^ exp 



il 



TkJcJ 



(2.14) 



1 /2 

with \h = {j§) being the magnetic length. 



M = Mo{9)\l + 2Re ^ A> [Af]* exp 



iAp ■ Ci (t) 



n 



(2.13) 



Having investigated the effect of the returning paths on 
the scattering process, we are ready to write an expres- 
sion for the dissipative current. Substituting Eqs. (2.13) 
and (2.14) into Eq. (2.8), we find 



Itt) 



coll 



de 



X M + 2Re ^ of, [af] * exp 



^£^4^ (a,C(t) • Ap) 

iAp ■ Ci (t) in{r - l)p]r\% i7r(r - l)e 



+ ■ 



+ 




(2.15) 



Equation (2.15) is the main qualitative result of this where 
section^^. It shows, that due to the presence of the 
returning orbits on one hand, and the largeness of the 
period on the other hand, the scattering process is ex- 
tremely sensitive to external fields applied to the system. 
As the result, a rich variety of the effects arises. The 
effects can be separated into two groups: (i) sensitive 
to the distribution function; and (ii) non-sensitive to the 
electron distribution. We will discuss those groups sepa- 
rately in following two subsections. 



Effects dependent on the form of distribution 
function 



Retaining only the contributions from the amplitudes 
with different winding numbers in Eq. (2.15), and consid- 
ering only linear terms in ^, we obtain familiar Shubnikov 
- de Haas oscillations Posc of the resistivity Pxx{B): 



.{B,T) 



Pxx{B = 0) 



, cos 



(2.16) 



. df{e) 
as — - — exp 



de 



i27rme 



huJr 



(2.17) 



are determined 



The form- factors, r]m — 

by the impurities at which scattering may occur during 
the circular motion of the electron, see Sec. II C for a 
more detailed discussion. The accurate expression for 
the factors ijm are written in Sec. IV B. 

Neglecting the effect of the applied electric field on the 
electron distribution function f{e), i.e. f{e) = /f(£) see 
Eq. (2.3), we obtain 



27r2mT 



(2.18) 



Therefore, the Shubnikov - de Haas oscillations arc ex- 
ponentially suppressed at high temperature T > loc/t^^- 
The consideration of effects non-linear in the applied 
electric field is more complicated task. Indeed, electric 
field may significantly change the distribution function 
/(e). In particular, the shape of the distribution function 
in a strong electric field may have nothing to do with the 
Fermi distribution. Because these effects are extremely 
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sensitive to the form of distribution function, their de- 
scription requires specifying the microscopic mechanism 
of the energy relaxation. 

We notice however that the smooth part of the elec- 
tron distribution function with the characteristic width 
Teff ^ T much larger than Uc results in exponen- 
tially small contribution to the resistivity. Therefore 
we focus our discussion to the high temperature limit 
T ^ fiLOc and neglect the Shubnikov-de Haas contribu- 
tion to the non-linear effects in further consideration. 
Unfortunately, the latter restriction does not allow us 
to avoid consideration of non-linear effects of the elec- 
tric field on the electron distribution function. Indeed, 
non-equilibrium component of the electron distribution 
function produced by the impurity scattering oscillates 
with period hujc-"^^ Substitution of such oscillating func- 
tion into Eq. (2.17) results in a contribution which is not 
exponentially small. An estimate^^ shows that this con- 



Therefore, all of the non- linear phenomena come from the 
effect of the external fields on the scattering cross-section 
itself. 

We consider the system under the effect of the dc elec- 
tric field only. In this case 

27r/ 

Ci{t) = dtC; (2.20) 

In the linear response |Ap(^;| ^ h, we obtain that the 
scattering rate is enhanced due to multiple returns: 

^^ = l + 2Re| „>[.<]•. (2.1) 

Postponing the estimate of the form-factors a^[af]* > 
until subsection II C, see also Sec. IV, we notice that 
Eq. (2.21) describes positive magnetoresistance. Indeed, 
it is intuitively clear that stronger the magnetic field, the 
larger the probability for electron to return to the same 
impurity. Thus, the contribution of the sum in Eq. (2.21) 
increases. 

One can see that the contribution of the periodic or- 
bit produces non- linear current voltage characteristics. 
Indeed, with the increase of the electric field \Ap<^i\/h 
becomes of the order of unity. Contribution of the re- 
turns with corresponding winding number I would be 
suppressed; electron after I turns simply misses impu- 
rity. If field is such |Ap^;^i|/7i ~ 1, then the contribu- 
tion of all returns will be suppressed. As the result, at 
dc fields larger than some value Eq the current voltage 



tribution may dominate all other effects, discussed below. 
A more detailed analysis of the non-linear effects on the 
electron distribution function and electron transport will 
be presented elsewhere.^'* 



B. Effects independent from the form of 
distribution function 

The contribution which survive the thermal or energy 
averaging are only those with the same winding number 
I = V . Retaining only such terms in Eq. (2.15), we can 
see that the energy dependence of the scattering cross- 
section vanishes, so that the energy integral can be eval- 
uated. The result does not depend on the distribution 
function anymore: 



(2.19) 

I 

characteristics becomes linear but with the slope deter- 
mined by the transport time in the absence of magnetic 
field, see Eq. (2.10). We estimate the value of Eq by 
noticing that according to Eq. (2.2), \dtC,\ = eEo/muJc, 
and \C,i\ — 2TTEo/meUj'^. Then, using |Ap| ~ 7i/^, where 
5 is the correlation radius of the impurity potential, we 
obtain the characteristic field 

eEo ^ mul^ (2.22) 

Accurate theory of nonlinear effects in the dc field is con- 
tained in Sec. V. 

Assume now that the ac microwave field with the fre- 
quency u) is applied together with the dc field. The ve- 
locity of electrons due to those fields dtC{t) and the dis- 
placement during I periods C;(t) = (^{t — 2itI/uJc) — C,{t) 
can be found as 

dtC{t) = Vdc + wCac cos Ojt, 

Clit) = 2C„eSin — cos - — . 

Wc \ ulc / \ Wc y 

(2.23) 

In the linear response regime those two velocities give two 
independent contributions to the current, i.e. dc response 
is not affected by ac radiation at all. Presence of the non- 
linear term in the collision probability Eq. (2.19) results 
in the photovoltaic effect. Indeed, expanding the expo- 
nent in Eq. (2.19) up to the second order and averaging 



dR 

dt 



coll 




dO eMo {6) (dtCit) ■Ap)ll + 2Re J2 [af ] * exp 
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the result over time, we estimate 



dtCit) exp 
dtCit) 
Vd 1 - 



n 



l-^(Ap.C.Wf 



Ap^C:c sin^ ^ 2'Kluj Ap^C^c sin ^ 



(2.24) 



It is important to emphasize, that at the frequency 
of ac field commensurate with the cyclotron frequency, 
u) = jiVc, the ac field does not affect the dc resistivity. 
This result is easy to imdcrstand by noticing that under 
this condition the impurity returns to its initial position 
during one cyclotron period. 

The last term in Eq. (2.24) represents the photo- 
voltaic effect and deserves some attention. Substituting 
Eq. (2.24) in Eq. (2.19) and neglecting the second term 
in Eq. (2.24), we obtain 



PxxjB) 
Pxx{0) 



^ / 2KlTrlljJ 2-kIw\ ^ r ^n* 

l + 2ReV 1 sin ] [af] , 

1=1 ^ 

(2.25) 

where m oc ICacP- The physical meaning of the photo- 
voltaic term is the rectification of the ac current due to 
the non-linear term in the collision integral. In the ab- 
sence of the dc field there is no preferred direction, so 
that the rectified current vanishes, whereas application 
of the dc voltage defines this direction. One can see that, 
due to the large factor 2ttuj/lUc 3> 1, this term can exceed 
unity even if |Ap(^^| < h. Therefore, the sign of the con- 
tribution from the returning orbits may be changed in 
comparison with the dc result, compare Eqs. (2.21) and 
(2.25). The zero- voltage dc resistivity may become nega- 
tive. Accurate results for the dc response in the presence 
of the microwave are collected in Sec. VI. 

Last comment concerns microscopic justification of the 
main assumption of Ref. 12, that even if the zero Re- 
current resistivity under microwave radiation is negative, 
it becomes positive at large enough applied dc current. 
Indeed, according to the arguments before Eq. (2.22), the 
contribution of the cyclotron orbits vanishes if the elec- 
tric field is large enough. On the other hand, these are 
the only contributions affected by the microwave radia- 
tion. Thus, at applied dc electric field exceeding Eq the 
current voltage characteristics becomes linear, with the 
slope determined by the transport time in the absence of 
magnetic field Eq. (2.10), and not sensitive to the effect 
of the microwave. 



C. Form-factors, self-consistent Born 
approximation, and classical memory effects 

Equation (2.19) derived on quite general grounds 
would describe all of the physics quantitatively if the 



finding of the form-factors ' were a trivial task. Un- 
fortunately, it is not so, and the quantitative analysis of 
the transport requires a machinery developed in the sub- 
sequent sections. The purpose of the present subsection 
is to explain the structure of the form factors qualita- 
tively, and to clarify the physical meaning and the status 
of the approximations to be employed. 

Let us start from the contribution of the returning 
paths with one turn, af '^. Those factors are disorder 
specific quantities and we have to average them over 
the disorder realizations. The estimate for such aver- 
aging proceeds as following. Let us assume that the 
two trajectories travel through the different impurities 
as shown in Fig. 3a. Then, every impurity scattering 
will randomize the sign of afa^. Therefore, the only re- 
maining contribution originates from trajectories which 
were not affected by other impurities during the cyclotron 
motion at all. The amplitudes can be estimated from 

|a<'>|2 = Po(f^), where Po{t) = exp{-t/Tg) is the 

probability for electron not to be scattered at an angle 
during time t, and 



1 



PFVF j d9M{9) 



h Ttr Ttr 



(2.26) 



is the quantum lifetime. This gives the estimate (which 
is actually an exact answer) 



af = tti = e 



(2.27) 



One may naively try to estimate the contribution from 
the paths with two turns from 



.<12 



which leads to 03=02= e~"^ However, this esti- 
mate is not correct. The reason is that the path which 
was scattered off an impurity on the first turn can be 
rescattered off the same impurity once again and give 
the contribution which is not random, see Fig. 3a. In 
the absence of the external electric field or the microwave 
illumination it gives 



1 - 



27r 



(2.28) 



which differs from the naively expected value. Moreover, 
it is clear that a2 will be affected by the motion of the 
impurity, i.e. the formfactors by itself are functions of the 
dc and microwave field [this effect was discussed from a 
different point of view in Ref. 29] . 

One can see, e.g. from Fig. 4b, that the larger the 
number of turns, the larger the number of the returning 
orbits must be taken into account. To perform analytical 
calculations, we utilize the self-consistent Born approxi- 
mation (SCBA) (first used for the disordered electrons in 
the magnetic field by Ando and Uemura^^). This approx- 
imation enables us to describe the combinatorial factors 



FIG. 3: a) Pair of typical semiclassical trajectories and; b) of 
those contributing to the classical memory effect. 



and the effect of the fields on the intermediate scattering 
processes correctly. 

The drawback of the SCBA is that the disorder poten- 
tials acting on electron on different semiclassical paths 
are assmned to be uncorrelated with each other. The 
justification of such an approximation requires that the 
typical distance between trajectories, see Fig. 3a, to be 
larger than the correlation radius of the potential ^. On 
the other hand this distance can be estimated as RcO ~ 
TiRc/{pfO- We obtain the inequality ^ ^ ^Rc/{pf£.) or 



hRc 

PF 



X2 



(2.29) 



where Xh is the magnetic length. Condition of the va- 
lidity of the self-consistent Born approximation was first 
established by Raikh and Shahbazyan^^ by an explicit 
calculation of the first correction to the self-consistent 
Born approximation. 

The case of weak magnetic field such that only single 
turn trajectories remain l/rtr <C Wc <C 1/Tg, requires 
additional consideration, even if the criteria (2.29) is sat- 
isfied. The trajectories close to each other, shown in 
Fig. 3b, give the main contribution to the interference 
of return trajectories, even though the fraction of these 
non-typical trajectories is small. Indeed, these trajec- 
tories travel through the same disorder. The scattering 
off the same disorder does not randomize the sign of the 
product Re(a^a{'*), and thus exponential estimate no 
longer holds. On the contrary, for 6 < one can 

roughly estimate from Fig. 3b 



af{e)a^{0) ~ exp 



27r / Ar 



where Ar is the typical distance between trajectories. 
Estimating |Ar| ~ RcSi, one finds for the angular inter- 
val of non-typical trajectories 



Rc 



\PfO \R 



„2 c2 \ 1/2 



\pf£.J \R. 



(uJcTtr) 



1/2 



The contribution of the scattering process to the trans- 
port scattering time (2.10), is proportional to 6f and the 




FIG. 4: Trajectories with (a) two and (b) three turns con- 
tributing to the form-factors aj^g; Inset: the semiclassical 
path of the electron in the vicinity of the impurity. 



typical scattering angle is fi/ippS,)- Therefore the rel- 
ative contribution to the non-typical trajectories to the 
change in all the dissipative processes, say resistivity, is 



Apa;a: 



C \ I n3/2 
) (WcTt,.) 



(2.30) 



This power law dependence should replace exp(— ^j^) 
dependence in all the formulas obtained in the self- 
consistent Born approximation for hu^^ S> T (not con- 
taining Shubnikov - de Haas oscillatory term). 

Finally, we notice that the product of two amplitudes 
for the electron propagation in the same disorder can 
be described as the classical probability of the circular 
path, and all the discussion can be recast into the notion 
of the classical memory effect (CME). It is not acciden- 
tal that the estimate (2.30) coincides with calculation of 
the CME magnetoresistance'^*' up to numerical prefac- 
tor. Accurate calculation of the non-linear effect within 
CME model will not be done in the present paper. It is 
important to emphasize, however, that the fact that the 
self-consistent Born approximation has to be corrected 
by classical memory effect, affects only overall prefac- 
tors in the non-linear effects and do not change the basic 
structure of Eq. (2.19).^^ 



III. QUANTUM BOLTZMANN EQUATION 
(QBE) 

In this section, we will derive the quantum Boltzmann 
equation. The propose of this derivation is to separate 
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the contributions remaining at 3> 1 from the very be- 
ginning. We use the standard Keldysh formahsm for the 
non-equiUbrium system^^. 



[idt - h) 



(1> 



A. Derivation of the semiclassical transport 
equation 

The derivation is very similar to that for the Eilen- 
berger equation^'*'^^. The matrix Green functions and 
the corresponding self-energies have the form 



G 



qR qK 





K 







(3.1) 



K 



where the component of the matrices are linear operators 

in time space and the one electron Hilbcrt space. 
The equation for the Green function is 

{idt -h)g= 1 + tG; G (idt - ^) = (3.2) 

where H is the one-electron Hamiltonian of the clean sys- 
tem, 1 is the short hand notation 't = Ik ® leS {ti — t2) 
and lK,e are the unit matrices in the Keldysh space and 
the one electron Hilbert space respectively. [Wc put 
?i = 1 in all of the intermediate formulas.] For retarded 
Green function we use 

{idt -h)g^ = S{t - ti)ie + t^G^; G^ {t < h) = 0. 

(3.3) 



and G^ = 



. For the Keldysh Green function it 

is convenient to take the non-diagonal component of the 
difference of two Eqs. (3.2) 



idt^H]:G 



(3.4) 

where [•; •] stands for the commutator. The next standard 
step is to separate the time evolution of the occupation 
numbers / and the wave-function of the system 



G 



K 



G^'-G' 



G^f - fG^ 



(3.5) 



In general / is an operator in both one-electron space and 
the time space. In the thermal equilibrium, however, one 
has simply 



fF = 

fF{e) 



ds 
2^' 



-ie{ti-t2) 



fFie) 



iT 



2sinh7rr(ti - ^2+^0)' 



(3.6) 



where T is the temperature in the energy units. Substi- 
tuting Eq. (3.5) into Eq. (3.4), one obtains the kinetic 
equation 

(dt + iH);f]=St{f}, (3.7a) 



= I + 



c 



(i) 



FIG. 5: (a)-(b) Self-consistent Born approximation (SCBA). 
Diagrams (c) are most important contributions not included 
in the SCBA, corresponding to the logarithmically divergent 
second loop weak localization correction^®, see e.g. Ref. 37. 
Diagram (c (i)) with one of the diffusons (d) reduced to one 
impurity line describes the first contribution to the classical 
memory magnetoresistance^" . 



with the collision integral given by 



st{/}= 



(3.7b) 



Next step is to write down a self-energy for the electron 
subjected to the random potential U (r) characterized by 
the correlation function 



(«7(ri)C/(r2)) 



/ 



W{q)e 



iq{ri-r2) 



For the sake of concreteness we will adopt the model with 

W{q) = W(0)e-«« (3.8) 

where ^ is the disorder correlation length. Equation (3.8) 
is an adequate description for the potential created by re- 
mote donors situated on the distance ^/2 from the plane 
of two-dimensional electron gas. Self-consistent Born ap- 
proximation, involving the summation of all the diagrams 
with non-intersecting impurity lines, see Fig. 5 is 



S = 



_fq_ 



W{q) 



(3.9) 



Self-consistent Born approximation is justified if two 
conditions 



^< Ah, 



27rn 



(3.10a) 
(3.10b) 



hold. Hereinafter Xh = {jg) is the magnetic length, 
and B is the applied magnetic field. The physical mean- 
ing of the condition (3.10a) was discussed in Sec. II, sec 
Eq. (2.29). Gondition (3.10b) allows us to neglect the 
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localization effects shown on Fig. 5c. Note, that the ex- 
ternal microwave radiation further suppresses the local- 
ization correction^*. 

We also assume small angle scattering 



(3.11) 



where is the Fermi momentum. This condition is not 
really essential for the physical processes but it allows 
for some technical simplifications.'*''' Moreover, based on 
Shubnikov de-Haas data, we believe that this regime is 
the most relevant for the experimental situation. 

We will consider the system in the classically strong 
magnetic field, so that the Hall angle u>cTtr is large (rtr 
is the transport time, see Eq. (3.50) below). The effects 
we will be studying are proportional to the inverse Hall 
angle and vanish in clean systems. Therefore, it is con- 
venient to solve time dependent problem for the clean 
system first, and then to consider the effect of the dis- 
order on the top of this solution. This program is easily 
accomplished by using the transformation Eqs. (2.1) and 
(2.2). Transformation (2.1) removes electric field from 
the Hamiltonian. Instead, the disorder potential becomes 
time dependent. We rewrite Eq. (3.9) as 



Cl2 = C(tl)-C(t2). 



'9Cl2 



(3.12) 



Separating the guiding center coordinate and the cy- 
clotron motion, we write 



P 



H = - n; r = R+ X^ip, 



(3.13) 



where n is the chemical potential and the operators obey 
the following commutation relations 



R,p 



0. 



The Green functions G{p;R) and the self-energies 
I](p; R) can be obviously written as functions of the op- 
erators (3.14). Using commutation relations (3.14), we 
obtain from Eq. (3.12) 



:Wi2 



^{p-'^)-J (2.)^ 
Wi2 = W{q)e-'1<^^ ; R{q) =R + X%iq. 

(3.15) 

Here we suppressed time indices. 

Next step is to separate the motion in the phase space 
into components parallel and perpendicular to the Fermi 
surface. For this purpose, we parameterize the cyclotron 
motion operators as 

1/2 

e^'^ + h.c; 
. 1/2 (3-16) 



1 



N ■ 



e'*' + h.c, 



where the integer 



U>r 



is introduced for convenience. To preserve the commuta- 
tion relations (3.14) for operators (3.16), the commuta- 
tion relation 



[n;0\ 



(3.17) 



is imposed. It follows from Eqs. (3.17) and (3.16) that the 
integer eigenvalues n > —Noi the operator h have the 
meaning of the Landau level indices. The Hamiltonian 
(3.13) acquires the form 

H = uje[n + S{fi)] ; d{fi) = l+N--f^. (3.18) 

Z (jJc 

All the previous manipulations were valid for any mag- 
netic field. Now we are going to make use of the large 
filling factor 



» 1. 



(3.19a) 



We will assume that the characteristic value of n con- 
tributing to the transport quantities is such that ]| n ||<C 
N, i.e. all the relevant dynamics occurs in the vicinity 
of the Fermi level. This assumption is justified provided 
that two conditions 



(3.19b) 



are satisfied, with Tq being the quantum elastic scattering 
time, see below. Those assumptions allow for the semi- 
classical consideration of the self-energy (3.12), which is 
presented below. 

Using Eqs. (3.19), we expand Eq. (3.16) as 



Px = Pf cos ip ■ 

Py = PF sin if - 



where 



1 

2i?; 
2^ [' 



/27V 



ne"^ + h.c. +.. 



PF 



X 



me 



Rc 



h.c] 



(3.20) 



/2iVA 



H 



H 



are the Fermi momentum and the cyclotron radius re- 
spectively. Substituting Eq. (3.20) into Eq. (3.15) and 
keeping in mind condition (3.11), we find 



t {h;^;R) = J '^^^ \U\\U^G{h;^;R{q)) t^lf/J 



C/j_ = exp iqRc sin (pq 



iq^RcSm2ipq q^R^ 



4pf 



Pi 



Un = exp ■ 



\fie'^Vq _ li e. 



ipq = (f-(pq, 



(3.21) 
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where q = q{cofi(j)q]sva(j)q). Matrix U±_ commutes with 
ip, i.e. in semiclassical sense it describes the scattering of 
the electron perpendicular to the Fermi surface. On the 
other hand, matrix t/y changes (p, i.e. it has the semi- 
classical meaning of the evolution parallel to the Fermi 
surface. Due to the condition (3.19a) those processes 
originate from parametrically different values of q and 
that is why they can be considered separately. We use 
the parameterization 



de 



Commutation relation (3.17) gives 



e-*''"G (6*, (^) , 



(3.22) 



Thus, using definitions of Eq. (3.21) we obtain 
ij^G (n; ^; R{q)) (}{ = j^J ^e-'^"G {6, <f) 

X exp •|i(7-Rc [siii(^g + 9) — sin<^g] 

iq^Rc[s.m{2ipq + 26) - sin2<^g] 1 
ApF I 



(3.23) 



The characteristic values of q entering into the integral 
can be estimated as qRc ^ -Rc/C ^ Rc/^h = V2N 1. 
Let us call the argmnent of the exponent f3{6,(j>). Be- 
cause, qRc ^ 1, the integrals will be determined by the 
saddle point determined by /3^(0, </>) = 0, and P'^{0, 0) = 
the latter condition gives ^ = 0. Thus, we write 

ei/3(M)«27^5(^^)<5[/3^(0,0)], 

in a sense that the saddle point integration on the LHS 
gives the same result as the integration in the RHS. Em- 
ploying this approximation in Eq. (3.23) and taking into 

account || n ||<C qRc, 



qRc 



S cos (fiq — 



q cos 2(fi^ 

2PF 



(3.24) 

where (5— function has to be understood in an operator 
sense. 

We substitute Eq. (3.24) into Eq. (3.21) and use 
Eq. (3.22) to find G{e = 0). With the help of the commu- 
tation relation (3.17) and using the small angle scattering 
condition (3.11), we obtain 



2ttV} 



^VF(g) exp [-igCi2 • 



X g 



<fi-\ ;-R- >^h1H 



cos I ifi 



q 

2PF 



sm I (p 



q 

2PF 



(3.25) 



where we introduced the analog of the Green function in 
the Eilenberger equation 

g (0; = ioJc^G (n + k; (p\ R^ . (3.26) 



Notice, that operator (p commutes with all other oper- 
ators entering into Eq. (3.25), and, therefore, it can be 
treated as a c-number. 

Equation (3.25) also can be rewritten in a different 
form 



St,t' = -ilCt,t'9{t,t';ip,R), 
where the kernel Kt t' is 



with 



W{q) 
A-k'^vf 



exp 



qd^ 

2pF, 



qd<p 

2pF, 



(3.27) 

dq, 

(3.28) 
(3.29) 



This form is more convenient for the further expansion 
in 1/^pF, as we show later. 

Closing this subsection we write down the expressions 
for the variation of the electron density, SN^, and the 
current density, j . Representing the coordinate operators 
by Eq. (3.13) and approximating the operators with the 
help of Eq. (3.20) we obtain 



SNe (r, t) = -me J ^^9^ [t, t; Vg] ; 



(3.30a) 



Rcsi{(p) - C{t), 



and 5g^ denotes the deviation of the Keldysh Green 
function from its equilibrium value in the absence of 
the external fields. The coordinate (^{t) is defined in 
Eq. (2.2). 

For the full electric current we have 



j{r,t) = eNc{r,t)dtC + j'-''Hr,t), 



(3.30b) 



where the first term in the current is dissipationless. At 
constant electric field it is nothing but the Hall current. 
The second term is given by 



/d 
—i{(p)g^[t,t;(p;rg]. (3.30c) 



The numerical coefficients in Eqs. (3.30a) and (3.30c) are 
written with account of the spin degeneracy, and A^e is 
the total electron density. 



B. Equation for the spectrum 

In this subsection we solve equation (3.3) with the 
Hamiltonian (3.18) and the semiclassical self-energy 
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(3.25) 



{idt -iOc[n + 5(m)]} G« {t, tv,n) = ^JLJll 



(3.31) 



Hereinafter, wc suppress ip and R arguments in the Green 
function and the sclf-oncrgy whenever they are the same 
in the both sides of the equations. Our purpose is to 
represent the in terms of the Green functions (3.26) 
only. 

For the calculation of the spectrum it is sufficient to 
keep the terms only to the zeroth order in small param- 
eter q/pF- Equations (3.25) and (3.27) then simplify to 



hi 



(Ci2g - »A|^ Vji) • ijip) 



(p;R 



(3.32) 



Here 



is nothing but the standard Born approximation for the 

quantum scattering time for small angle scatterers, and 
the dimensionless function 



h\ {x) = 



JdqWiq) 



1 + 



(3.34) 



characterizes the effect of the external electric field dur- 
ing the cyclotron motion between electron returns to the 
same impurity. 

We will look for the self-energy in the form 



iS^(i,ii) = 



2to 



+ 



1=1 



(3.35) 

where the coherence factor A describes the phase accu- 
mulation during one period and it is defined as 



A = exp 



z27r(5(/i) 



c ' q 



(3.36) 



and iS/* are to be found self-consistently. The time shift 
operator is defined as 



TH = t 



2nl 



(3.37) 



We look for the solution of Eq. (3.31) in the form 



iG^ {t,ti;n,ip) 



-iwci5((ii) (t-ti ) g~ "27^ 



27r 



xe{Th-t^e(ti-f^+^t^ |. 



(3.38) 



Substituting Eqs. (3.38) and (3.35) into Eq. (3.31), and 



using e 



i2'Kln 



1, e*^"9Je 



ip + x we obtain the chain 



of equations 

(t,fi;^) = <?f (ti;^) 

- f dt2Y, [^""'*2; 'fi + UJc {t2 - tl)] Gtm it2, tv, <p) 

(3.39) 

where 



gPit,ti) = gp_,(ft,ti) 



(3.40) 



with the initial condition — 1. 

Final step is the self-consistency procedure which 
amounts into substitution of Eq. (3.38) into Eq. (3.26). 
It gives^^ : 



Ht-h) 



■oc- 

Y^X^gP{ti)s(f't-ti) 



(3.41) 



where A is defined in Eq. (3.36). Using Eq. (3.41) in 
Eq. (3.32) and extracting Si of Eq. (3.35), we find 



S^{t,^) = hi 



{Ci{t)e-iXj,VR)-i{^) 



9^ (t, f) 



We use the short hand notation 

c,it) ^ c {f-H) - at) , 



(3.42) 



(3.43) 



where the finite time shift operator is defined in 
Eq. (3.37). 

Equations (3.39) - (3.42) constitute a complete system 
for the spectrum averaged over disorder. Note, that the 
Green function Gi is expressed in terms of Qm with m < I 
only. 



C. Equation for the distribution function 

In this subsection we will reduce Eqs. (3.7) to the 
canonical Boltzmann form. According to Eq. (3.25) the 
self-energies T, do not depend on h. It suggests that 
the distribution function / docs not depend on n either. 
This observation enables us to substitute Eq. (3.5) into 
Eq. (3.26) and perform the summation over k with the 

help of Eq. (3.41) and relations 5^ = - [g^]\ f = p. 
We find 
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^ OO 

2 [9'' -9^- 9"] = / {t, t';^,R) + Y^ A'^f {f't; ip, r) f (f% t' , ^, r) 

1=1 

00 

+ E (^0* / ^. r) [9? (^'*'; ^> 



(3.44) 



where the coherence factor A is defined in Eq. (3.36), and the time shift operator is given by Eq. (3.37). 

Substitution of Eq. (3.44) into Eqs. (3.30a) and (3.30c) yields the connection between the distribution function / 
and the observables: 



5N, (r, *) = 2me y ^ I 5f {t, t; ip,rg) + 2Re A'^f (fh; if, r^) / (f't, t; cp, r^) I ; (3.45a) 



t) = 2epF / ^^(^) I / (*' rg) + 2Re ^ A'gf (f% >p, Vg) f (f% t; <p, rg) | , (3.45b) 



where we used the short hand notation Vg = r — Rcei{(p) — <^{t). 

To derive the kinetic equation, we substitute / into Eqs. (3.7). Equation (3.7a) gives 



d_ d_ 
dt^df 



d_ 

dip 



f{t,t';ip,R)=St{f}t,t 



(3.46a) 



According to Eq. (3.7b), the coUision integral is defined in terms of the electron self-energy, the latter is given by 
Eq. (3.27). Substituting Eq. (3.27) into Eq. (3.7b) and using the relation Eq. (3.44), we obtain the following expression 
for the collision integral: 

St{f}t,t'=-{lC{t,t')fit,t')] 

^A'l []C{t,fh)gf' (fh)] f (f%t') - [/C(i,i')/ {f%t') {fh)] | 

^(A*)' / {t,fH') [t{fH',t')gf (fh')] - [t{t,t')f (t,fh') gf (fh') 



+ 



(3.46b) 



1=1 



In Eq. (3.46b) kernel fC is given by Eq. (3.28), functions 
gf''^ are defined in the previous subsection by Eq. (4.1), 
the coherence factor A is defined in Eq. (3.36), and the 
time shift operator is defined in Eq. (3.37). Note that 
we suppressed ip and R arguments in the entries of 
Eqs. (3.46b) for brevity. The first Une in Eq. (3.46b) 
corresponds to the classical scattering off an impurity; 
the second and the third lines describe the retarded in- 
terference corrections due to the returning orbits. 

One property of the kinetic equation (3.46) is worth 
emphasizing because it is a crucial check of the consis- 
tency of the approximation we made. Consider the con- 
stant electric field E, so that 



m = ~ct—. (3.47) 



Then the distribution function 

f{t,t') = fF{t-t')e-''^'^^'-''^ 

where the equilibrium distribution function is given by 
Eq. (3.6), null both the collision integral and the left 
hand side of Eq. (3.46a). In the energy representa- 
tion /(e) = fpis + eER) corresponds to the thermody- 
namic equilibrium of the system in the moving coordinate 
frame. 

For the small angle scattering, see Eq. (3.11), we ex- 
pand the rotation operator exp{qdip/2pf) in the powers 
oiq/pF and obtain the following expression for the kernel 
Eq. (3.28): 

t = + ICj + IC^ + ko (3.48) 
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where the first term contains only even angular harmon- 



ics: 



in Eq. (3.49b) to display explicitly that the operator Kj 
is Hermitian as guaranteed by the relation 



Ttr 



X hi 



(3.49a) 



dhz{x) 
dx 



2ixh2{x). 



the second term contains odd angular harmonics: 



Ttr 



Ttr 



(3.49b) 

and the third term represents the angular diffusion: 



K.^ = dtph2 

Ttr 



d^. (3.49c) 



The remaining term ko describes contributions which are 
of the order of l/{pF^)^ smaller. In Eqs. (3.49) we intro- 
duced the transport mean free time 



1 1 f dq 11 wiun some eiiecLive e/j 

— = - 5" / —q^W(q)o=i — - — — , (3.50) teraction prevents Te f f 

Ttr ^-KVfPIJ 27r Tq{pF£,r Tn tlnP nf il 



Finally, the first term iC±_ of the kernel Eq. (3.48), see 
Eq. (3.49a), is responsible for the evolution of the dis- 
tribution function perpendicular to the Fermi surface, 
which does not mix angular harmonics of the distribu- 
tion function with different parity. Similarly to KLj, in 
the first line of the collision integral Eq. (3.46b), K,± 
describes the classical effect of the electric field on the 
electron distribution - Joule heating. The other terms in 
Eq. (3.46b) describe the returns of the cyclotron trajec- 
tories, which may result in oscillating components of the 
distribution function with period TilUc-'^'^''^^ We remark 
that term Eq. (3.49a) of the collision integral can not be 
considered separately from inelastic processes such as the 
electron-electron and electron-phonon interactions. In- 
deed, taking the zeroth angular harmonics of Eqs. (3.46a) 
and (3.49a), one finds the correction to the distribution 
function, that grows infinitely in time. This is just a sig- 
nature of the energy absorbed by the system from the ex- 
ternal field. Elastic impurities alone can not stabilize the 
distribution function. The electron-electron interaction 
suppresses large deviations from the Fermi distribution 
with some effective Te// whereas the electron-phonon in- 



from an infinite increase. 



and the dimensionless function 



/l2 {x) = 



Jdqq^Wjq) 
Jdqq^W{q) 



l-3x^ 



{l + x- 



\3 ■ 



(3.51) 



Function hi{x) is defined by Eq. (3.34), and i = 
(cos ip, sin (p) . 

Let us discuss in more details the meaning of the com- 
ponents of the kernel Eq. (3.48). The third term in 
Eq. (3.48), given by Eq. (3.49c), describes the angular 
diffusion. Its contribution to the collision integral sup- 
presses angular harmonics of the distribution function 
rather than the zeroth harmonics. 

The second term ICj in Eq. (3.48), see Eq. (3.49b), 
represents the scattering process accompanied by simul- 
taneous creation of odd-angular harmonics and energy 
shift. It is this term that is responsible for the dissipa- 
tive current. When its is substituted into the collision 
integral Eq. (3.46b), the first line of the collision integral 
describes an instantaneous scattering and gives the clas- 
sical conductivity. The second and third lines describe 
the interference due to the returns of the cyclotron tra- 
jectories. All non-linear effects considered in the follow- 
ing sections originate from the second and third lines of 
the collision integral Eq. (3.46b) with the full kernel IC 
replaced by ICj. The linear in term in Eq. (3.49b) 
may be safely omitted since it does not give rise to any 
effects relevant for the future consideration. We kept it 



In the remainder of the paper we will consider only the 
phenomena associated with fCj, that are not sensitive to 
effects of the external field on the distribution function. 
Therefore we neglect the contribution to the collision in- 
tegral Eq. (3.46b) originating from the term Eq. (3.49a) 
of the kernel Eq. (3.48). This contribution may be ne- 
glected if the energy relaxation time is small, the condi- 
tion which is generally not valid. In recent preprint^"^ an 
estimate of the contribution to the dc resistivity from os- 
cillating component of the distribution function was pre- 
sented. When inelastic processes are weak, this contribu- 
tion is larger than the contributions studied in the present 
paper. Nevertheless, the effects considered here are ro- 
bust and are described by different system parameters, 
therefore they deserve a separate consideration. The con- 
tribution which depends on the form of the electron dis- 
tribution function will be presented elsewhere.^'* The im- 
plicit assumption everywhere will also be (T, Tg//) <C fi, 
i.e. the electron system is degenerate. 



IV. LINEAR TRANSPORT 

The purpose of this section is two-fold: (i) to demon- 
strate how the solution of the QBE is obtained for the 
simplest case and to reproduce relatively known results; 
(ii) to derive formulas for the spectrum which can be used 
as building blocks for consideration of more elaborate ef- 
fects in the further sections. 
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We begin with the solution of Eq. (3.39) for the spec- 
trum. In the hncar regime, we can put ft-i = 1 in 
Eq. (3.42). After this simphfications entries of Eqs. (3.39) 
become independent of the angle ^p and time ti . We find 
with the help of Eqs. (3.40) and (3.42) 



^1 m=l 



(4.1) 



9i 



« - gp (0) = g^^. 



with the initial conditions = 1, g^^ = 0. 

Non-linear recursion relations (4.1) can be resolved ex- 
actly with the result 



^ , , Ll[{x + l)a] ^ a{l - x)Ll, [{x + l)a\ 
yi{x) = r— ; 1- 



91 



l + l 

I 



X = 



27r' 



/ + 1 

277 



(4.2) 



where L["(.t) is the Lagucrre polynomials^. 

Next step is to find the distribution function. To do 
so we use Eq. (3.6) for equilibrium distribution and solve 
Eq. (3.46a) to the leading order in l/{uJcTtr) ^ 1 and in 
the first order in C(0: ^^e Eq. (3.47). In this approxima- 
tion only the collision term (3.49b) contributes. Taking 
the limits f2 ^ ii — ^ for integer I, we find 



nr \Pt + '^cJ I 



sinh 



27r^r(;-m) 



m=0, m^l 



sinh 



(4.3) 



where C,i{t) is defined in Eq. (3.43). 

Equations (4.2) and (4.3) are sufficient to calculate 
the linear response of the electric current Eq. (3.45b) 
to the applied electric field within the self consistent 
Born approximation at arbitrary temperatures and mag- 
netic fields. First we discuss the high temperature limit, 
T ^ TiLOc, when the conductance is a smooth function 
of the applied magnetic field. Then we take into ac- 
count l/u)c oscillations of the conductivity which appear 
at T/uc^ 1 (Shubnikov-de Haas oscillations). 



A. Ac- and do- transport at high temperatures 

At T » only the first term in Eq. (4.3) remains and 

all other terms arc exponentially suppressed. Introducing 
Ne =Pp/{2tt) and substituting Eqs. (4.2) and (4.3) into 



Eq. (3.45b), we find with the help of Eq. (3.36) 



{df +Lj'^y \meTtr 



27r/ 



{ft) 



(iA) 



where the time finite shift operator is given by Eq. (3.37). 
First term in (4.4) describes usual scattering contribution 
and the subsequent terms result from the multiple returns 
to the same impurity. In the frequency representation 
Eq. (4.4) may be written as 

j(f) = a(^)(^)i;., 

where a'^'^'^Lo) is the conductivity tensor: 



{—iuje + LOc) f e'^Nf 



nieTtr 



1 + 2^. 



2Trl 



(4.5) 



At = 0, Eq. (4.5) gives the non-oscillating correction 
to the diagonal dc resistivity ^xx- Using LVcTtr ^ 1, one 
writes p^x = crxx{<^ = 0) [pxy] and Eq. (4.5) yields 



where 



Pxx{B) = porjo 



PD 



2tt 



Ttr 



is the Drude resistivity, and 

oo 

Vo{a) = l + 2j2^[Lli{a)Y 



-al 



The asymptotic behavior of Eq. (4.8) is 



%(a) 



e-'"(l-a)' +0{ 



-3a\ 



at a 1 (weak magnetic field), and 

16 fl 
37r V a 



(4.6) 



(4.7) 



(4.8) 



(4.9) 



(4.10) 



at a <C 1 (strong magnetic field). Function 770(0;) to- 
gether with its asymptotes is plotted in Fig. 6. 

Equation (4.6) was written to the leading order in 
l/{uJcTtr)- To obtain the correction to the Hall coeffi- 
cient one has to solve Eq. (3.46a) and take into account 
the term (3.49c) in the first order perturbation theory. 
The Hall coefficient can be expressed in terms of the 
third powers of the coefficients (4.2). The final result. 
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FIG. 6: The solid line represents the high temperature mag- 
netoresistance curve Eq. (4.6), a = 2tx j (lOcTq). The dashed 
curve in the main panel is the high field approximation (4.10). 
The dashed curve in the inset is the low-field asymptotic ex- 
pression (4.9). 



At weak field, a ^ 1, function Ti{w, a) is well described 
by the first few terms: 

Ti {w,a) w l-|-2e""cos27rw-|-2e"^"(l-Q!)^cos47rit;-|-. . . . 

(4.15) 

At strong magnetic field, a ^ 1, we use the asymptotic 
expression of the Laguerre polynomials in terms of the 
Bessel functions Jn{y)^^ and obtain 



-lx/2 



I 



-L\-i (l^) 



Ji (2/Vi) 



1 



(4.16) 



Substituting Eq. (4.16) in Eq. (4.15) and employing the 
Poisson summation formula we obtain for a <^ 1 



Ti {w, a) 
where 
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E ^1 



A:— — oo 



7r\w — k\ 



(4.17) 



37:9 {2 ~ \x\) 



dy cos xy 



Mx) 



(2 + a;)6l(2- |a;|) 



however, will have the smallness l/(wcTtr)^ and that is 
why we will not write down the explicit form of those 
corrections. The Hall coefficient in this approximation is 



Pxy 



B 
ecN^ 



1 



,2^2 



(4.11) 



where !Fh{x) vanishes exponentially at a; — > 0. 

At finite frequency w, Eq. (4.4) gives, in particular, the 
oscillatory dependence of the absorption of microwave 
radiation with field E{t) — ReE^e"'^'^^ on frequency w. 
We find 



Re El,a^^{uj)E^ = 

,2 



to:. 



rrieTtr 

^ — 2lljuJc cos 2/3 w 27r 



(c.2-^2)2 



(4.12) 



where parameter (3 describes the polarization of the field 
by parameterization of E,^ as 



V2E^ 



e+ cos /3 -I- e_ sin /3, e± = e ± iee, 

(4.13) 



where e is the unit vector. For the circular polarization of 
the microwave /3 — 0,±7r/2. For the linear polarization 
along e, (3 — tt/A. 

The dimensionless function J-'i('w,a) = J-'i{w + l,a) 
represents the normalized coefficient of microwave ab- 
sorption: 



-la 



C0S2'kIw n -I ^n2 ,x 

— [Ll^ (la)] . (4.14) 



X <^ (i + x^) E 



-AxK 



Functions K{x) and E{x) are complete elliptic integrals 
of the first and second kind respectively, and function 
TLi{x) can also be obtained as a convolution of two semi- 
circle densities of states. Equation (4.12) and its asymp- 
totes (4.15) and (4.17) are consistent at \uj — uJc\Ttr ^ 1 
with the result of recent preprint. 




W 



FIG. 7: Dependence of normalized microwave absorption co- 
efficient J^i{w,a), see Eq. (4.14), on frequency w = uj/uic at 
three values of a = The function J^i{w,a) is periodic 

in w with the period 1. Note, that at strong magnetic field 
(a = 0.3) function vanishes at intervals around w = k + 1/2 
with integer k. 

According to Eq. (4.17), at sufficiently strong mag- 
netic field the absorption coefficient J-'i(w,a) vanishes 
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at frequency intervals such that |u; — A;| > l-^/aj-K with 
integer fc, as one may expect for the case, when the den- 
sity of states has gaps between Landau levels, see e.g.^''. 
Numerical investigation oi T\(w^oi) at intermediate val- 
ues of a ~ 1 allows us to find the threshold value of 
magnetic field, when the gap appears in the two level 
correlation function within the SCBA; this value of the 
magnetic field corresponds to a « 0.65. Figure 7 shows 
T\{w,o) for three values of a, including the threshold 
value a ~ 0.65. We note that the vanishing of the two 
level correlation function at some energy interval is an 
artifact of the SCBA. However, the correction to this re- 
sult (tails in the density of states) are known*^ to drop 
exponentially with the increase of the Landau level index 
N and we disregard such tales in our study. 




B. dc transport at low temperature 

At low temperature terms in the second and third lines 
of Eq. (4.3) become important. Substitution of these 
terms into Eq. (3.45b) yields for the resistivity at uJcTtr ^ 
1, compare with Eq. (4.6), 

Pxx{B,T) 

PD 

(4.18) 



where po is the Drude resistivity, Eq. (4.7), and 



Tix) 



sinh(27r^a;) 



, T(0) = 1. 



The disorder dependent coefficients r]i{a) — ri-i{a) are 
given by 

^ / «(|/ + fc| + |fc|) ^ 

m (a) = 2^ exp I 1 [a)g\i+k\ (a), 

k=-oo ^ ^ 

(4.19) 

with gi (a) defined in terms of the Laguerre polynomials 
by Eq. (4.2). 

Term with ^ = in Eq. (4.18) reproduces the smooth 
part of the magnetoresistance (4.6) and / > 1 represent 
the Shubnikov-de Haas oscillations. The asymptotic be- 
havior of function 77; (a) is the following. At low fields, 
a ^ 1 only the first few terms are relevant 

r,i = 2e-"/2 -)- 2(1 - a)e-3"/2 + Q (e^^^/sj ; 
7^2 = (3 - 2a) e-" + O ; 

= (4 - 8a + 3a2) e'^^/^ + O (e-^^/^) ^^'^^^ 
m=0 (a'-ie-'"/2) . 



FIG. 8: The harmonics -qi of Shubnikov-de Haas oscillations, 
Eq. (4.18) as functions of magnetic field and the quantum 
scattering time, a — -22—. 



At high magnetic field a <C 1, we use Eq. (4.16) and 
obtain for l^/a <C 1 



77; (a) 
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1 - 



2l\ 



0{l 



4 2\ 

a 



(4.21) 



Coefficients ??/(a) obtained from Eqs. (4.19) and Eq. (4.2) 
are plotted in Fig. 8. 

Below in this paper we assume that temperature is suf- 
ficiently high, T ^ LOc- This assumption allows us to dis- 
regard the Shubnikov-de Haas oscillations in transport 
quantities. 

In this Section, we applied the quantum Boltzmann 
equation to calculate the linear response of electron sys- 
tem on the applied electric field. The approach developed 
here enabled us to describe the resistivity at arbitrary 
values of the parameter u>cTq. Our findings are in accord 
with the results of Ref. 1,27,40 and differ by an overall 
numerical factor from the corresponding result of Ref. 41. 
Particularly, the strong magnetic field asymptote for the 
smooth part of the magnetoresistance matches the result 
of Ref. 27. The amplitude of the Shubnikov-de Haas 
oscillations calculated in this Section is consistent with 
the previous analysis of magnetoresistance oscillations in 
Refs. 40,41. We also derived an expression for the absorp- 
tion rate of microwave radiation, Eq. (4.12). Asymptotes 
of our expression in weak and strong magnetic fields coin- 
cide with the results presented in Ref. 23. Having made 
sure that the consequences of the quantum Boltzmann 
equation are consistent with the results obtained by dif- 
ferent methods, we will proceed with the description of 
electron transport beyond the linear response. 
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V. NON-LINEAR dc EFFECTS 

A strong dc electric field produces non-linear effects 
on (i) the even harmonics of the distribution function, 
see Eq. (3.49a); and (ii) the elastic scattering processes 
(spectrum), see non-linear terms in Eq. (3.49b). The 
first mechanism roughly corresponds to the heating and 
it strongly affects system properties determined directly 
by the electron distribution function, such as Shubnikov- 
de Haas oscillations of the conductivity, see Sec. IV B. 
As we have already noticed in Sec. IIIC, the distribu- 
tion function is very sensitive to the details of inelastic 
processes. If, however, the temperature is large. 



(5.1) 



from the very beginning wo can restrict our considera- 
tion to the non-linear effects on the electron scattering 
process. In the remainder of the paper we consider only 
the high temperature limit Eq. (5.1). 

Similarly to Eq. (4.3), we take into account only colli- 
sion term (3.49b). Neglecting exponentially small terms 
and assuming UcTtr '> 1, we solve Eq. (3.46a) and obtain 



2wSf {T%t;<fi) = 



V 

IVcTtr 



C{ti,t2;ip) 



dt2 



Ci2 • 



(5.2) 



We notice, that the distribution fmiction 5f{TH,, t: (p) 
does not depend on time t. In the linear response regime 
/i2 = 1, see Eq. (3.51), and Eq. (5.2) reduces to the first 
term in Eq. (4.3). 

Using Eq. (3.45b) and Ng = Pp/^n we find 



2eN, 



271- ^ 




+ 



OO 

2 J2 I Ap'sf {fh; ^) d-' [fff {fh; ^) C {TH, t; ^) 



(5.3) 



Equation (5.3) is more complicated than its linear re- 
sponse counterpart, because the spectrum, determined 
by gf^it, if), depends on the applied dc field. It prevents 
one from using Eq. (4.2); therefore Eq. (3.39) should be 
solved again. Using Eq. (3.42), and Eq. (2.2) for the 
dfE = 0, we obtain instead of Eq. (4.1) 

e«(t,(p)=gf(^) 



Eq 



(5.4) 



where hi is defined in Eq. (3.34), and we introduced a 
scale for electric field 



Eq 



2-Ke 



(5.5) 



Explicit angular dependence in Eq. (5.4) makes the so- 
lution for an arbitrary magnetic field difficult. We con- 
sider only limiting cases of weak, uJcTq <C 1, and strong, 
WcTq ^ 1, magnetic fields. 

At weak magnetic field, ujcTq <C 1 we can keep only first 
two non-trivial terms in Eq. (5.3). Solutions of Eq. (5.4) 
for Z = 1, 2 are angular independent: 



31 = 50 = 1; 52 = 1 - 



El 



27r 



1/2 



Substituting these functions into Eq. (5.3) we obtain the 
solution in the form 



, E 27r 

^2 TT, 



\E\ \rn,,,Ujf.Ttr J " ^ \ Eo' UJcTq 



(5.6) 



where the dimensionless function ^2(2;, a) in the weak 
magnetic field can be expanded as 

^2(a;,Q!> 1) = a;|l-|- 

2 (1 - 8x2) g-2. 



+ 



(1 + Ax'^f' 




For the weak dc fields, \E\ <C Eq, Eq. (5.7) match 
Eq. (4.9). 

At strong magnetic field WcTq » 1 the second angular 
harmonics in the solution of Eq. (5.4) is suppressed by 
a factor of a = 27r/(a;cTq) in comparison with the zero 
angular harmonics. Neglecting this correction and intro- 
ducing 

~9^{E/Eq, a) = gf'iE/Eo, a)e-"'/2 (5.3) 
we obtain from Eq. (5.4) 



a 
2J 



Equation (5.3) simplifies to 



91 {x,a)-a ^ . 
„=i {l + w?x^)' 

(5.9) 



^ (1-2ZV) , 



tr(l + /^x2)= 

(5.10) 

find a solution of Eq. (5.9) as a correction to Eq. (4.2) 
and use approximation similar to Eq. (4.16): 



9i{x,0i) 



y^Ji(2V^) El {lx;af); 



1, 2/<i; 

V'J^ y ^ 1 



(5.11) 
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Substituting this expression into Eq. (5.10), we find with 
the logarithmic accuracy for x 1 
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y 



4a 



In 



f- 



(5.12) 



The second term in Eq. (5.12) represents the suppression 
of the renormahzed transport time due to the apphed 
electric field, compare to Eq. (4.10). 

At strong electric field, \E\ > ^\i/2 , Eq. (5.9) can 



be solved by perturbation theory in l/{LOcXq) 

(I — m) 



g^{x, a) ^ 1 



'i (1 + im?x^) 



1/2 



(5.13) 



This gives with the help of Eq. (5.10) and Poisson sum- 
mation formula 



T2{x,a) 



27ra 

^2 



Sir' 
.t3/2' 



\/a<x<l. (5.14) 



For the strongest fields \E\ ^ Eq the main contribution 
to the current becomes linear in field: 



T2ix,a) 



,C(3)-2aC(2) 



a; > 1, 



(5.15) 



where (^{x) is the ^-function. 

Dependence of ^^2(2;, a) on x is plotted in Fig. 9 for 
several values of a. Function J-2{x,a) is calculated ac- 
cording to Eq. (5.10) with functions gi{x,a) obtained 
from recursive equation Eq. (5.9). At strong magnetic 
field, !F2{x,a) exhibits a non-monotonic behavior with 
a minimum &t E ^ Eq. At strong electric field, E 3> 
Eq, all curves approach the zero magnetic field result, 
!F2{x,a) = X, since the strong electric field destroys the 
interference effect of electron motion along cyclotron or- 
bits. 

Figure 10 shows the asymptotes of function T2{x,a) 
for three intervals of the strength of electric field, see 
Eqs. (5.12), (5.14) and (5.15). For comparison, we also 
show T2{x, a) calculated directly from Eq. (5.10) for a = 
1/40. 




FIG. 9: Nonlinear dependence of the dissipative current on 
the applied electric field x = E / Eo at high-temperature 
TijjJc <^ T, for diflterent values of the magnetic fields, a — 

2-E juJcTq. 




FIG. 10: The plot shows three asymptotes of function 
^"2(2;, a) for a = 1/40: i) for small a; < 1 [Eq. (5.12), dashed 
line]; ii) for -^/q <C x [Eq. (5.14), solid line]; iii) for large 
a; ^ 1 [Eq. (5.15), dot-dashed line]. The dotted line repre- 
sents function T-iix^a) at a = 1/40 calculated directly from 
Eq. (5.10). 



VI. EFFECT OF MICROWAVE RADIATION ON 
dc TRANSPORT 

Consider the two-dimensional electron gas in a mag- 
netic field LOcTtr 3> 1, subjected to a monochromatic mi- 
crowave radiation together with the dc field 

E{t) =E + Re£;„e-*"*, (6.1) 

where E^^ is a complex vector in the plane of two- 
dimensional electron gas. 

For the strong magnetic field such as the filling factor 
is small, 1/ ~ 1, the effect of microwave radiation was 
considered in detail in Ref. 9, and the linear response 



for the short range disordered was analyzed in Ref. 11. 
Our goal is to extend these studies to the small angle 
impurity scattering and to the non-linear dc response. 
The first direction will make the theory more adequate 
for the description of the experiments^^^'^, whereas the 
second development provides the microscopic grounds of 
the theory of zero-resistance state^2_ Tj^g latter issue is 
analyzed in further details in the subsequent section. 

To characterize the microwave power in dimensionless 
units, we introduce 

j^^ ^c (t^g+^'-2a;^eCOs2/3) /cj.y E^ ■ E*^ 

(6.2) 
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where characteristic field, Eq, is defined in Eq. (5.5). 

The polarization of the microwave is described by the 
angle /3 and the unit vector e as prescribed by Eq. (4.13). 
We will introduce also the parameter 



, , / CO + U!c „ 

'yicoj = arctan tan p 



(6.3) 



which describes an elliptic trajectory of a classical elec- 
tron in the magnetic and microwave fields. 

The only difference of Eq. (6.4) from Eq. (5.2) is the 
time dependence of the distribution function and the 
spectrum due to the oscillating microwave field. Work- 
ing in the approximation of large Hall angle ujcTtr ^ 1 
and large temperature T ^ hojc, we once again solve 
Eq. (3.46a) with the collision term (3.49b). We obtain 



{dt + ujcdip) 



C {ti,t2;ip) 



d_ 



i{'p)C{ti,t2)h2 



where C(^i;*2) = C(^i) ^ C(^2) with vector defined in 
Eq. (2.2), and A is introduced in Eq. (3.36). Explicitly, 



C(ii,i2) _ uj,{t2 ~ h)eE_ ^ ^^^^{t2 - ii) 



2TTEn 



X Re e (e+ cos 7 + e_ sm 7) e 2 



The spectrum of the system depends both on the mi- 
crowave radiation and the applied electric field. From 
Eqs. (3.39), (3.40), and (3.42) we find similarly to 
Eq. (5.4) 



9l 



\t) = g^{t,t) = gti {f-h,t) 



t-tl 



{t, h,ip)= gf^{t, ip)--J2 / dt2g^{t2, f + t2U>c)GfLrn (*2, t, if) hi {Zi,m{t2, h)) , 



Zl,m{t2,tl) =i{ip + t2U)c) 



'1 m=l-^0 

mE 



(6.5) 



-pr- + VP sin Re (e+ C0S7 + e_ sin7) e-'^^^i+'^'e -c 

Eq IjJc 



where hi is defined in Eq. (3.34), time shift operator is given by Eq. (3.37), field Eq is defined in Eq. (5.5), dimensionless 
power of the microwave radiation V is given by Eq. (6.2), and the angle 7 is given by Eq. (6.3). 

We substitute the electron distribution function Eq. (6.4) into Eq. (3.45b) and obtain the following expression for 
the dc component of the electric current in terms of gf-{t,ip), compare to Eq. (5.3): 



UlcTtr 







(6.6) 



and (. . . )t stand for the time averaging over the period 

of the microwave field. 



A. Weak magnetic field, oJcTq <C 1 



Equation (6.6) together with the recursion relation 
Eq. (6.5) determines the electric current to all orders both 
in the microwave power V and dc electric field E/iJo. 
We remind that our consideration is valid for large filling 
factors V \, large Hall angle ojcTtr ^ 1, large tem- 
peratures T ^ hujci and under the conditions of the ap- 
plicability of self-consistent Born approximation (3.10). 
Further simplifications are possible for certain limiting 
cases, which will be considered below. 



In this case we can limit ourselves with only first non- 
trivial term in Eq. (6.6). Because gi = 1, further cal- 
culation is reduced to straightforward angular and time 
integration in Eq. (6.6). The regimes where the compact 
analytic results are available are listed below. 
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1. Circular polarized microwave radiation 



For the linear response in dc electric field E, we find 



nieOJ^Ttr 



(6.7) 



where dimensionless microwave power is defined in 
Eq. (6.2) and 



J^3{'P,w,a) 



e-" (2-Psi 

1 + 

(l + Psin^ 



sm TTW 



.5/2 



-Ve "sin27rw- 



(4 - -P sin^ nw) 
1 +7'sin2 ^^y/^ 



+ 0(6-2"). 
(6.8) 
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-■- 75 = 0.5 
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w 



At 75 = 0, Eq. 



matches Eq. (4.9). 



Structure of Eq. (6.8) deserves some additional discus- 
sion. The second term in brackets describes the effect of 
the microwave on the elastic scattering process. In the re- 
gion of the applicability of the theory uJcTq ^ 1 , its value 
can never become larger than the first term and their sum 
is always positive. This follows from the fact that the 
elastic transport cross-section is positive by construction 
no matter what kind of rcnormalization it acquires. The 
third term is the photovoltaic effect discussed in Sec. II. 
Its sign depends on the frequency of the radiation to and, 
remarkably, on the power of the microwave radiation. It 
is noteworthy, that this term may make the current flow 
opposite to the electric field even at small magnetic field 
due to the presence of possibly large factor Qttuj/uJc in 
front. Finally, we emphasize non-monotonic dependence 
of the photovoltaic effect on the microwave power. The 
frequency dependence of the dc resistivity at weak field 
is plotted in Fig. 11. The corrections to the Hall coeffi- 
cient are small as l/{ujcTtr)'^ and will net considered here 
explicitly. 

If the microwave power is small 



^ . 9 

VSIT? < 1, 



FIG. 11: Frequency (w — ui/ujc) dependence of jF3(7-',w,a) 
at fixed value of the microwave power parameter V from 
Eq. (6.2) [notice, that it corresponds to the actual microwave 
power dependent on the frequency ui] in the weak magnetic 
field uJcTq — O.Stt. The upper panel presents curves corre- 
sponding to the weak microwave field, "P ^ 1, whereas the 
solid curve in the lower panel corresponds to the strong mi- 
crowave radiation 7^ = 4. Curves are calculated for the circu- 
larly polarized microwave. 



where 



- SVe^" sin'' nw 



2 (1 - 2x^) e' 
1 H 

(1 + x^f^ 



„ . 2 4~27x^ + 7x^ 



SiTwVe "sin27rti;- 



2(l + x2)^/2 
x^ 



2(l + a;2)^/2 
+ OCp2e-°) + 0(6-2") 



(6.10) 



one can expand Eq. (6.8) up to the first order in V. In 
this case, the whole non-linear dc response affected by 
the microwave can be found, compare with Eq. (5.6) 

\E\meUJ^Ttr \Eq UJ^ UJcTq J 



Function !F^{x,V ,w,a) is plotted in Fig. 12. One can 
see, that at large electric field x ^ 1 the Ohm law is 
restored and the microwave radiation becomes irrelevant 
in accord with the conjecture of Ref. 12. We will discuss 
consequences of negative values of JF4 in Sec. VII. 

At weak dc field, \E\ ^ £'0 one can keep only terms 
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X 



FIG. 12: Non-linear dependence of J^4{x,'P,w,a) on the 
strength x — \E\/Eo of the dc field for different values of 
power parameter V and for u = 7.25tJc. Compare with Figs. 9 
and 11. Curves are plotted for the regime of weak magnetic 
field ujcTq = TV and for circularly polarized microwave. 



linear in dc field in Eq. (6.9), 



nieUJ^Ttr 



6Ve~ 




sm 



l + 2e" 



TiLO sm ■ 



(6.11) 



It coincides with the corresponding expansion in 
Eq. (6.8). 



2. Arbitrary polarization of the microwave radiation 

The compact results can be obtained for the first or- 
der expansion in Vsiu^— < 1. Polarization of the 
microwave is characterized by the parameter 7 from 
Eq. (6.3). We find 




FIG. 13: Dependence on microwave frequency of the electric 
current for a linear polarization of the microwave field: i) 
e||i5 (dashed line) and ii) e_LiS (solid line). The curves are 
calculated for weak magnetic field Ucfq = tt and weak power 
of microwave V = 0.2. The strength of the electric field is 
\E\ — O.5E0. The plot shows that the possibility to form a 
current in the direction opposite to the electric field depends 
on the polarization of microwave field. 



where functions for the dipole and quadruple angular har- 
monics are given by 



'iiTwsa\2'KW 3(1 — 42:^) sin^ TTCJ 

(6.14a) 

TTw sin 27rz« (3 — 4x^) sin^ ttw 
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[l + x"^)"^ 



[l + x^f 



(6.14b) 



for P < 1, a > 1. 

For the current linear in the dc field and bilinear in the 
microwave field, Eq. (6.12) simplifies to 



j''^(7)=jt^+<5jaSin27, 



(6.12) 



where j^'*' represents the isotropic component of the 
current and coincides with the current produced by 
circularly polarized microwave field, Eq. (6.9). The 
anisotropic component is given by 



meOJ^Ttr 



E -2 [eE) 



Eo 



J 1 T) 

5 I 77^ ' ' ' 
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(6.15) 



We emphasize that the anisotropy of the electric current 
versus the applied electric field appears both in the linear 
and non- linear dc transport. 

In Fig. 13 we plot the dc resistivity for the linear po- 
larization, /? = 7r/4, of microwave field for the cases E\\e 
and EJ-B. One can see from Fig. 13 that the condition 
for the electric current to flow against the applied elec- 
tric field E depends on the polarization of the microwave 
field with respect to the direction of E. 
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B. Strong magnetic field, ujcTq ^ 1 

In this case, Eq. (6.5) can also be significantly simpli- 
fied. We will limit ourselves with the first order expansion 
in microwave power V. 

First we analyze the first order correction in V to the 
spectral functions in Eq. (6.5). By inspection, one 
can see that this equation contains terms either slowly 
changing during the cyclotron period or oscillating with 
frequencies ncu. The oscillating with frequency 2uj term 
do not contribute at all to the final answer, whereas the 
term oscillating with frequency uj can be taken into ac- 
count perturbatively. Thus, for the redefined spectral 
functions according Eq. (5.8), we obtain analogously to 
Eq. (5.9) 
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where the angle independent component 5/^ satisfies the 
following recursion relation 



FIG. 14: Nonlinear dependence of Ti, Eq. (6.22), on the 
strength x — \E\/Eo of the dc electric field at several val- 
ues of the microwave frequency w = uj/uic in strong magnetic 
field LOcTg = 407r at P = 0.05. 
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2^ Eo~gtm9^ J, , A 

(6.17) 

and the angle dependent component 5g^^{'p) can be 
found from 

(6.18) 

The initial conditions for the recursion relations arc 
gQ = 1, 6gQ — 0. Above we introduced the short-hand 
notations 



Si(m) = sin^ 



2El - m^E^ - 3to2 sin 27 [2{eEf - E^] 



A{E^ + m^Ey 



(6.19) 



and 



„ , , . TTmuj . f d(j) 
'^2(™) = sm Ef\ / — exp 



2TTiuJ(f> 



" (6.20) 



[mE ■ [(e+ C0S7 + e_ sin 7) • 2(0)] 



E^ + (mi(0) • E)' 



where 7 is given by Eq. (6.3). 

One can see, that due to the oscillating factor e 
in Eq. (6.18) the contribution of this term is suppressed 
by additional factor of lOc/^ in comparison with contri- 
bution of Si in Eq. (6.17). Nevertheless, even Si, which 



describes the effect of the microwave radiation on the 
density of states, is suppressed in comparison with pho- 
tovoltaic effect by a factor of uJc/^- Thus, in the consid- 
eration of the transport at 



CO ^ LJc 



(6.21) 



we replace gf {t,Lp), defined by Eq. (6.17), with gl^{t,if ) 
obtained from Eq. (5.9). All the further formulas of this 
Section are valid in this high-frequency limit only. 

We present the dc current in the form of Eq. (6.12), 
where the polarization dependence is characterized by 
factor 7, Eq. (6.3). Keeping in mind condition (6.21), we 
obtain from Eqs. (6.4) and (6.6) the following expressions 
for the functions defined in Eqs. (6.9) and (6.13) 
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(6.22) 



and 
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J^r,(x,V,w,a) = SttwV } ^ zn::KH9i [x,a}; 



Tl{x,V,w,a) 
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2 tii^ + i^x^) 



7/2 



(6.23a) 
(6.23b) 



forP<Cl, ?i'^l,Q!<Cl Here, the spectral functions gi 
are solutions of Eq. (5.9). 

In Fig. 14 we present J-'i(x,V,w,a) as a function of 
the strength of the electric field x — \E\/Eq for several 
values of frequency w = to/uc at strong magnetic field. 
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where for I a; I < 2 



1 + ^ ^Z^^2 



2 



(6.26) 



and 7^2(2;) = otherwise, see also discussion in the last 
paragraph of Sec. IV A. Function 7^2(2;) has the mini- 
mum at Xrriin ~ 0.834, where Ti.2{xmin) ~ 0.726. Cor- 
respondingly, function J-'Q{V,w,a) becomes negative if 
the microwave power V exceeds Vmm ~ 0.29v^/w (here 
a 1 and w ^ I). This expression demonstrates that at 
strong magnetic fields, a ^ 1, already a weak microwave 
is sufficient to create a state with zero-bias negative re- 
sistance. 



FIG. 15: Dependence of J-%{V ,w,a) on frequency w = uo/u)c 
for several values of V. The upper panel represents the curves 
calculated for a — 0.65, and the lower panel shows the curves 
for a = 0.30; cf. Fig. 7. 



i^cTq = 407r. We observe that the effect of microwave field 
on dc current is significant only in the non-linear region, 
\E\ < Eq. At stronger electric fields '> Eq, the effect 
of microwave on the dc current disappears, see Sec. II for 
a discussion. 

We also calculate the dc current linear in the dc field 
E and bilinear in the microwave field. For simplicity 
we consider only the isotropic component of the current, 
which survives at 7 = in Eq. (6.12) and corresponds 
to the current produced by the circular polarization. In 
this case we use the spectral function gi{x,a), given by 
Eq. (4.2) and obtain 



-Te{'P,w,a), 



(6.24) 



where 



Tg{V, w, a) = lim 



■ =%(a)H — 2~aw 

(6.25) 

and functions ?7o(a) and T\(w^a) are defined by 
Eqs. (4.8) and (4.14) respectively. Relation between the 
absorption spectrum and the microwave frequency de- 
pendence of the photovoltaic effect (6.25) was argued re- 
cently in Ref. 18 on the basis of a "toy model" . Depen- 
dence of J-q{V a) on frequency w — uj/lUc is shown 
in Fig. 15 for several values of V. [Note that fixed V 
corresponds to the frequency dependence of the actual 
microwave power, see Eq. (6.2).] 

At strong magnetic field we use Eqs. (4.10) and (4.17) 
to find the asymptotic form of function J-'e(V,w,a) at 



VII. FORMATION OF INHOMOGENEOUS 
PHASES AND CURRENT IN DOMAINS 

Results of the previous section qualitatively consistent 
with the conclusions of Ref. 9-11 indicate that there is a 
region in the parameter space where the linear dissipative 
conductivity becomes negative. According to Ref. 12, 
spatially homogeneous state of such system is unstable 
and break itself into the domains characterized by zero 
dissipative resistivity and conductivity and by the clas- 
sical Hall resistivity, see Fig. 16. In the analysis of such 
state one can ask two main questions: (i) what the spatial 
structure of the domain wall and the boundary conditions 
fixing the position and the size of the domains are; (ii) 
what the values of the current and the electric field inside 
domains are; the value of electric field can be found by 
the local probe measurement. 

First question has to be answered by analyzing spa- 
tially inhomogeneous problem by taking into account the 
gradient term in Eq. (3.49b) and the Poisson equation; 
this question is left for future study. Here, we use the 
results of Sec. VI to address the second question.^* 

To clarify further consideration, let us discuss the rela- 
tion between applied current and voltage in more details. 
In all of the above analysis we assumed that the electric 
field E is applied and the current j is measured, the cur- 
rent has both the dissipative and Hall components; the 
corrections to the Hall coefficient are small as 1/{uj'^t^^). 
Restoring the Hall current we write the expression for 
the total dc current up to the terms ©(cj^^r^^^) 



e^N.Eo 



^, E ^ LU 271 \ E 1.^ 

T{—,V,—, — eE, 

Eq UJc UJcTqJ \E\ p^y 



(7.1a) 

where p^y is the classical Hall coefficient, see Eq. (4.11), 
and J- is the tensor defined for different situations in 
Eqs. (5.6), (6.9), and (6.13). [Tensor structure appears 
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c) 



FIG. 16: Domain structure^^ for the Corbino (a) and Hall bar 
(b) geometries. For the Hall geometry, the applied current I < 
jdL is accommodated by the shift of the domain wall without 
any voltage drop, Vx — 0- At I > joL {V > Ed) for the Hall 
bar (Corbino) geometry domain structure is destroyed and 
the state with finite dissipation is stable (c). Note, that the 
construction (b) does not describe the current pattern near 
the leads. 



due to the microwave radiation with polarization other 
than circular, see Eqs. (6.12) - (6.15), (6.23). ] Mi- 
crowave power is characterized by Eq. (6.2) and the field 
Eq is given by Eq. (5.5). This relation is convenient to 
use for the Corbino disk measurement scheme. For the 
Hall bar geometry Eq. (7.1a) can be easily inverted 



E = -JoPdsT ( —,V,—, 

Jo i^c 



2tt \ ij 



-g-^+Pxyij, (7.1b) 



where = —1, pn is the classical Drude resistivity (4.7), 
and 
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Eq 

Pxy 



siV, 



27r 



(7.2) 



is the electric current scale for non-linear effects. Equa- 
tions (7.1) shows that both Hall bar and Corbino mea- 
surements should exhibit the similar non-linear proper- 
ties, as will be discussed below. 

We mainly focus our consideration on the circular po- 
larization of microwave radiation; non-circular polariza- 
tion is briefly discussed in the end of this section. Then, 
tensor from Eqs. (7.1) is reduced to scalar and the con- 
dition of the local stability of the state takes the form^-^ 



where x is the solution of 



(7.3a) 



, u) 27r 
0; d^T{x,V,—, )>0. 



(7.3b) 
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FIG. 17: Upper panel: Phase diagram of a 2DEG in weak 
magnetic field (a = 2) m V — w = ^ coordinates [1 = 
(y = 0) for the Hall bar (Corbino) geometries]. Region (A) 
is the dissipative state; region (B) is the zero resistance (con- 
ductance) state. Lower panel: Phase diagram for oj = 4.3a;c 
(dashed line on the upper panel) in 'P — J coordinates for the 
Hall bar {V — V for the Corbino) geometries. The zero re- 
sistance (conductance) state exists if the current through the 
Hall bar (voltage drop between edges of the Corbino disk) 
does not exceed the critical value Ic (Vc). The same curve 
defines the value of the spontaneous current jd ~ Ic/L (elec- 
tric field Ed = Vc/L) in domains. The line labelled by 
is the numerical solution of Eq. (7.6) and the line "2-3" is 
a schematic interpolation beyond the linear expansion in mi- 
crowave power. Points Vu.i are the same as in the upper 
panel. 



All the further analysis is reduced to the substituting of 
appropriate limit of the function (6.9) into the stability 
condition (7.3). 

The regime of the weak magnetic field a — ^ l is 
simplest. According to Eq. (6.7), the zero current state 
is stable if JF3 > 0, thus, the equation 



^3 



V,—,a 

UJr 



= 



(7.4) 



gives the boundary between dissipative and zero resis- 
tance state (ZRS) for the Hall bar geometry or zero 
conductance state (ZCS) for the Corbino disk geome- 
try iriV — ijJ plane, where ^P, a » 1^ is given by 

Eq. (6.8). The curve given by Eq. (7.4) is plotted in the 
upper panel of Fig. 17. For w = ^ ^ 1 the analytic 
estimates for the "phase boundary" fines are 



Vi 



dnw sin 2'kw 
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sm nw 



Vi « 1; (7.5a) 
5"/'Pi, Vu » 1. (7.5b) 
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The zero resistance state is impossible not only at too low 
microwave power but also at excessive microwave power 
(reentrance transition). Indeed, a weak microwave radi- 
ation does not produce strong enough photovoltaic cur- 
rent to compensate the dissipative current. On the other 
hand, as we discussed in Sec. II, strong microwave radia- 
tion suppresses the electron returns to the same impurity 
and thus destroys the non-linear effects. 

At microwave power within the zero-resistance region 
Eq. (7.3b) has the solution at a; ^ 

:f.UvA,^)=o. (7.6) 

For the low microwave power response J^4 is given by 
Eq. (6.9). 

Phase boundary given by Eq. (7.6) is shown on the 
lower panel of Fig. 17 by the 1 — 2 — 3 curve. In the 
vicinity of the lower boundary (segment 1 — 2) and at 
CO ^ Wc we have 

where Pi is given by Eq. (7.5a). As the power increases, 
the current in domains reaches maximum and then de- 
creases. This non-monotonic behavior is schematically 
shown by the 1 — 2 — 3 line in the lower panel of Fig. 17. 
The corresponding segment (2 — 3) may be obtained from 
calculations outlined in Sec. VI A 1 beyond the bilinear 
response in the microwave field, which was not done in 
the present paper. [Based on the results presented in 
Sec. VI A 1, only point 3 of this segment is known.] How- 
ever, there is no reason to expect any singular behavior 
of this curve. 

The lower panel of Fig. 17 may be also used as a phase 
diagram in V — I plane, where / = jL is the total 
current through the Hall bar of width L or in "P — V 
plane for the Corbino geometry, see Fig. 16. 

We now turn to the discussion of the strong magnetic 
field regime cUcTq '> 1. Naively, one would expect that 
increase of the magnetic field would change the phase 
diagram of Fig. 17 only quantitatively by rearranging the 
boundary line. However, this expectation is not correct. 
We start from the phase diagram on V^w plane. Fig. 18a. 
The condition for the boundary between the dissipative 
and ZRS (ZCS) (7.4) is modified as 

^5(^,^,0^=0, (7.8) 

where is given by Eqs. (6.25) and (6.26). Solution 
of Eq. (7.8) gives the line 1 - 2 - 3 in Fig. 18a. On 
c;an sec that the region of the instability shrinks with 
increasing of the magnetic field. It is not the end of 
the story though. According to Fig. 14, see curves for 
w = 5.15 and w = 5.25, the state with positive zero-field 
resistance but with dissipative electric field antiparallel 



to the electric current at some finite current is possible. 
The boundary line (curve 3 — 2 — 4 — 5 in Fig. 18a) for 
such state is given by, see also Eq. (7.3b) 



d,J^4(x.,V,—,^] =0, (7.9) 

where is given by Eq. (6.22). The solution of Eq. (7.9) 
is shown as the line 2 — 4 — 5 in Fig. 18a. Therefore, the 
phase diagram becomes more complicated. The region of 
the ZRS (ZCS) has the same properties as its counterpart 
for the weak field. On the other hand, in the coexistence 
region (C), see Fig. 18a, c, both the homogeneous dissipa- 
tive state with zero current and the domain structure of 
the ZRS (ZCS) are locally stable. We believe that such 
bistability can cause the hystcrctic behavior of the I —V 
characteristic of the sample, see Fig. 18(e). 

The aforementioned complication translates into the 
qualitative change in the phase diagrams in V—V (V—I) 
coordinates, see Fig. 18 (b,c) in comparison with the 
lower panel of Fig. 17. Equations for the lines on 
Fig. 18(b,c) are 

5x^4 (x,V, ^, ^) = 0; Lines V,,2 {ha) 
^4 {x,T, ^, ^) = 0; Line (h), 

and V = xLEq (J = xLjo). The physical meaning of 
those lines is illustrated on Fig. 18f. 

The region D in Fig. 18 (b,c) represents the state with 
negative differential conductance (for Corbino geometry). 
In this case, the homogeneous state is unstable, whereas 
the zero-resistance state is not possible. The instability 
of the homogeneous state, see Fig. 19 leads to the forma- 
tion of the domain structure with the charge distribution 
similar to the Gunn domain^^. This structure will be 
moving from the boundary to the boundary with the ve- 
locity $^(jJcT* /{2'k), thus the domain will be annihilated 
on the contact with the other one formed on the opposite 
contact; so that the current pattern will be oscillating in 
time rather than stationary. 

Finally, let us discuss the role of the anisotropy of the 
dissipative conductivity tensor in the formation of zero- 
resistance (conductance) state for the linear polarization 
of microwave. Here, two situations are possible, (i) both 
main components of the linear resistivity tensor are neg- 
ative though different; (ii)the main components of the 
linear resistivity tensor are of different signs, see Fig. 13. 
Study of the regime (i) can be reduced to the previously 
studied case by rescaling of the coordinate, currents and 
field, such that equations V-j=0, Vx£l = are kept 
intact. It docs not change the state of Ref. 12 qualita- 
tively, though extra singularities may be needed to ac- 
commodate the change in the boundary conditions. For 
case (ii), the homogeneous state can be shown to be un- 
stable, whereas the domain structure with closed current 
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FIG. 18: (a) Phase diagram of a 2DEG in strong magnetic field {a I) in P — w — coordinates [/ = (V = 0) for the Hall 
bar (Corbino) geometries]. The dissipative (A), the zero resistance (conductance) state (B) and the coexistence regions (C) are 
shown. (b,c) Phase diagrams for lu/ljc = wi and uj/lJc = W2 (vertical dashed lines on the upper plane) in V — I coordinates for 
the Hall bar (V—V for the Corbino) geometries. Lines Vz{Ii) also describe the current in domains as functions of the microwave 
power. The dissipative instability region is denoted by D. (d,e) The I — V , (V — I) characteristics for the Hall bar (Corbino 
disk) geometries a.t V = 1^2 and V = 1^4 respectively; (f) Relation of the position of the border lines to the results of Fig. 14. 




FIG. 19: The schematic picture of the current distribution in 
the instability region D for the Corbino geometry. 



loops would violate the condition § E-d\ — 0, because on 
such contour there must be regions of positive resistance. 
The details of current pattern for this case requires the 
further investigation, we believe, however, that the sta- 
tionary solution for this case is not possible and domains 
oscillating in time will be formed. 



VIII. CONCLUSIONS 

In this paper, we derived the kinetic equation within 
the self consistent Born approximation for large filling 
factors. The obtained equation are written in terms of 
the Green functions integrated in the phase space in the 
direction perpendicular to the Fermi surface similarly to 
the Eilenberger equation for normal metals and super- 
conductors. Our system of equations takes into account 
the effect of electric and magnetic fields on the elastic 
scattering process, i.e. on both the spectral function and 
the electron distribution function. 

Armed with the quantum kinetic equation for the limit 
of large Hall angle, we described the following phenom- 
ena: (i) dc and ac magnetoresistance in the linear re- 
sponse; (ii) non-linear dc current-voltage characteristic; 
(iii) influence of oscillating microwave electric field on 
dc current. It is important to emphasize that the non- 
trivial effects of the theory are described in terms of only 
two free parameters, time Tq which can be extracted for 
the Shubnikov-de Haas oscillations, and the characteris- 
tic electric field Eq from Eq. (5.5). The major problem of 
the presented paper is the lack of the consideration of the 
inelastic processes and consequently, effects related to the 
form of electron distribution function. The treatment of 
inelastic processes will be presented in ref. 24. 

We conclude by mentioning another consequence of the 
proposed in Ref. 12 domain mechanism of zero resistance 
(ZRS) and zero conductance states (ZRC). Namely, ac- 
cording to our finding, the zero dissipative current rep- 
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resent the interplay of two effects: elastic scattering off 
impurities and the photovoltaic effect; electric field in the 
domain is found from the condition that these two effects 
compensate each other on average. However, those pro- 
cesses are statistically independent. Consequently, this 
statistical independence of two processes may be revealed 
through the current noise in the ZRS or ZCS, which is 
not expected to have any singularity in this regime. The 
analysis of this noise can be performed by slight modi- 
fication of the equations derived in the present paper in 
the spirit of e.g. Ref. 43 and left as a subject for the 
future research. 
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it will not be important for the qualitative arguments, see 



Sec. Ill for the rigorous theory. 

Reliable estimate for the short range scatterer can be ob- 
tained by substitution ^ — > Ti/pf, Tq = Ttr in the final 
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The non-linear effects of the electric field on the electron 
distribution function, which are not taken into account, 
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ferent phases of the electron system. However, the overall 
structure of the phase diagrams remains intact even if the 
non-linear effects are considered. 



